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Summary 


Introduction 


A semiempirical method is presented for the esti- 
mation of aerodynamic forces and moments acting on 
a steadily rotating airplane. The airplane is divided 
into wing, body, and horizontal and vertical tail sur- 
faces. The strip theory is employed to wing, body, 
and horizontal tail surfaces to determine their con- 
tribution to the total aerodynamic coefficients. The 
strip theory takes sectional, two-dimensional static 
aerodynamic coefficients as input and predicts rotary 
aerodynamic coefficients. For wing and horizontal 
tail surfaces, such input data are usually available 
from static wind-tunnel tests. However, for the body, 
the static test data are not normally available in the 
required form. Therefore, a semiempirical approach 
is developed for the prediction of static aerodynamic 
characteristics at combined high angles of attack and 
sideslip to produce the required input to the body 
strip theory. For the vertical tail surface, strip theory 
is not applied because of its low aspect ratio; there- 
fore, a direct method is employed. All the mutual 
interferences between wing, body, and tail surfaces 
are ignored except the mutual interference between 
horizontal and vertical tail surfaces (shielding effect). 
A mathematical model for estimating the shielded 
area of the vertical tail at high angles of attack is de- 
veloped based on available experimental data. The 
predictions of the strip theory based on static aero- 
dynamic input data were found to differ considerably 
from the corresponding rotary balance test data. For 
wing and horizontal tail surfaces, the rotational flow 
model proposed by McCormick (NASA CR-165680) 
is used to determine additional contributions of these 
components. For the vertical tail, a model called 
“secondary flow effect” is developed to account for 
the curved flow effects over the vertical tail. 

The analysis of the paper is quite general. How- 
ever, some specific applications are considered in this 
paper. This theory is applied to a light, low-wing, 
single-engine, general aviation airplane, extensively 
studied in the NASA Langley stall/spin research pro- 
gram. 

In general, the predicted trends of the rotary 
(spin) aerodynamic coefficients are consistent with 
the spin-tunnel rotary balance test data. The agree- 
ment is reasonably good at low and moderate an- 
gles of attack and moderate spin rates. However, at 
higher angles of attack or high spin rates, the differ- 
ences become significant, particularly in estimation 
of side force and yawing moment. These deficien- 
cies are attributed to various difficulties in modeling 
complex aerodynamics of the spinning airplane. 


The prediction and analysis of airplane stall/spin 
characteristics have been of great interest to design- 
ers since the beginning of aviation. This problem 
has assumed more importance in recent years be- 
cause the uncontrolled motions of the airplane as- 
sociated with stall/spin have caused significant loss 
to the military and civil aviation (ref. 1). Although 
considerable progress has been made in recent years 
in the areas of experimental and flight testing tech- 
niques related to stall/spin problems, an adequate 
mathematical model is still not available for simula- 
tion of aerodynamic forces and moments in spin. 

The phenomenon of loss of damping in roll and 
autorotation of stalled wings has been generally rec- 
ognized as a primary cause of spin for light airplanes. 
The spinning motion usually involves sideslip, pitch- 
ing, rolling, and yawing about the body axes, and the 
airplane descends towards the Earth with its center 
of gravity describing a helical path. In steady-state 
spin, the aerodynamic forces and moments acting on 
the airplane are in equilibrium with inertia forces 
and moments. However, the establishment of bal- 
ance between aerodynamic and inertial moments has 
been found to be of crucial importance in attaining 
a steady-state spin. If the airplane cannot find such 
a balance, the spinning motion remains oscillatory. 
The radius of spin is usually of the order of 1 semi- 
span of the wings. However, with an increase in angle 
of attack a and spin rate, the spin radius approaches 
zero. In steep spin (a s t a ll < ot < 45°), the dominant 
part of the aerodynamic moment is the rolling mo- 
ment which mainly comes from the wings. In flat spin 
(45° < a < 90°), the major component of the aero- 
dynamic moment happens to be the yawing moment, 
and the chief contribution comes from the body as 
the vertical tail is often rendered ineffective because 
of the shielding effect created by the wake from the 
horizontal tail. Additional information on the air- 
plane spin can be found in references 2 and 3. 

The aerodynamics associated with a spinning air- 
plane is quite complex. It is dominated by extensive 
regions of separated flows and complex mutual inter- 
ference between various aerodynamic surfaces. Still, 
it remains as one of the most challenging problems 
of aerodynamics. The traditional approach based on 
the stability derivative concept does not prove useful 
in the analysis of spinning motion because the aero- 
dynamic forces and moments at poststall conditions 
display high nonlinearity with angle of attack and 
spin rate. This problem, at present, is still beyond 
the reach of computational fluid dynamic methods. 
In view of this, a semiempirical analysis based on as 
much physical information as possible and derived 



from experimental work can play an important role 
in the theoretical analysis of airplane spin. 

Glauert (ref. 3), Gates and Bryant (ref. 4), and 
Wykes, Casteel, and Collins (ref. 5) employed strip 
theory to estimate the aerodynamic characteristics of 
a steadily spinning airplane. The strip theory takes 
static aerodynamic data as input and predicts the 
rotary aerodynamic coefficients. These studies were 
limited to the analysis of steep spins. The problem 
of flat spins, where both the angle of attack and spin 
rates can be quite high, was not given much atten- 
tion. Nonavailability of an adequate aerodynamic 
mathematical model led to the development of the 
spin-tunnel rotary balance apparatus for generating 
the pertinent aerodynamic test data (ref. 6). By in- 
corporating the spin-tunnel rotary balance aerody- 
namic test data into the six-degree-of-freedom equa- 
tions of motion, it has been demonstrated that the 
steady-state spins can usually be predicted (refs. 7 
and 8). Thus, the experimental data base generated 
by rotary balance tests has played a very useful role 
in understanding the stall spin problems. On the 
other hand, for a spinning airplane, the experimen- 
tal investigations dealing with flow visualization and 
pressure measurements similar to those carried out 
for rotating propellers (ref. 9) and which could pro- 
vide insight to develop a comprehensive theory are 
almost nonexistent. 

In the present study, an effort has been made to 
estimate the aerodynamic forces and moments acting 
on a light spinning airplane. Here, the airplane is 
considered to be steadily rotating about a vertical 
axis analogous to a model undergoing spin-tunnel 
rotary balance tests so that the present theory ties 
in as closely as possible with the rotary balance test 
data. Also, this condition approximates quite well 
an airplane in a steady-state spinning flight. Angles 
of attack up to 90° and reduced spin rates up to 1.0 
are considered in the present analysis. 

The airplane model is considered to be divided 
into wing, body, horizontal and vertical tail surfaces. 
The effect of power is ignored. The strip theory 
used for wing and horizontal tail surfaces is based on 
the input of sectional, two-dimensional static wind- 
tunnel test data (ref. 10) and is similar to the meth- 
ods used in earlier studies (refs. 3 and 4). How- 
ever, for the body, the static wind-tunnel test data 
in sectional coefficient form are not usually avail- 
able. Therefore, a semiempirical approach is devel- 
oped to predict the static aerodynamic characteris- 
tics at combined high angles of attack and sideslip 
to generate the required input to body theory. All 
the interference effects between wing, body, and tail 
surfaces are ignored except the mutual interference 
between horizontal and vertical tail surfaces (shield- 


ing effect). A procedure for estimating the shielded 
area of the vertical tail is developed based on experi- 
mental data and is incorporated in the calculation of 
vertical tail aerodynamic characteristics. 

The predictions of the strip theory based on the 
static aerodynamic input data were found to differ 
considerably from the corresponding rotary balance 
test data (refs. 11 and 12). This discrepancy is at- 
tributed to the presence of curved or rotational fluid 
flow. Herein the term “rotational flow” is used to 
refer to the curved nature of the flow due to an- 
gular velocity fi. The simple strip theory accounts 
only for the kinematics associated with angular veloc- 
ity n and has no mechanism to handle the physical 
aspects of the flow pattern. A brief description of 
the flow field over rotating airplane surfaces is pre- 
sented based on available information in the litera- 
ture. For wing and horizontal tail surfaces, the rota- 
tional flow model proposed by McCormick (ref. 13) 
is employed to determine additional contributions of 
these components. The rotational flow effects are ig- 
nored for the body. For the vertical tail, a secondary 
flow model is proposed to account for its curved flow 
effects. 

The strip theory is applied to a light, single- 
engine, low wing general aviation airplane which has 
been extensively studied as model A in references 10, 
11, 12, 14, and 15. We have presented calculations 
for the basic configuration of the model A airplane 
of references 10, 11, 12, 14, and 15, whose tail con- 
figuration is designated as tail 4, and two other tail 
configurations (3 and 5). The presentation of results 
follows the pattern of spin-tunnel rotary balance test 
data of references 10 through 12. The nomencla- 
ture adopted is also similar to that used in these 
references. 

For application of the present method to other 
light airplane configurations, it is necessary to have 
the input empirical data in the form similar to that 
discussed. 

Some partial results of this investigation have 
been presented by the authors to recent conferences 
(refs. 16 and 17). An extension of the strip theory 
approach to wings of arbitrary planform is presented 
in reference 18. 

Symbols 

A axial force 

AR aspect ratio 

horizontal tail span 
local width of body cross section 

by vertical tail span 

&W wing span 


2 


C A 

axial-force coefficient, - T xia (, ^ orc ? 

jP^Sw 

Cd 

drag coefficient 

Cl 

lift coefficient, y - 77^ — 

^pUooSyy 

Ci 

rolling-moment coefficient, 
Rolling moment about eg 

7>pU£oS\\'b\v 

Cm 

pitching-moment coefficient, 
Pitching moment about eg 

Cn 

normal-force coefficient, 
Normal force 

^ pUoc^W 

Cn 

yawing-moment coefficient, 
Yawing moment about eg 

^pU^Syjbw 

Cx 

two-dimensional axial-force coefficient 

Cx,c 

axial-force coefficient of two- 
dimensional square section 

Cy 

side-force coefficient, f°J ce 

^pUooSw 

Cy 

side-force coefficient of unshielded 
vertical tail 

Cy,c( B) 

side-force coefficient of idealized two- 
dimensional body cross section of 
model A airplane 

CyA SC) 

side-force coefficient of two- 
dimensional square cross section 

c 

local chord 

c d 

drag coefficient of circular cylinder 

c f 

chord of flat plate 

c w 

wing mean geometric chord 

D 

drag force 

d f 

maximum width of ffat-plate wake 

Fy,v 

vertical tail side force 

h y 

vertical distance between X-axis and 
centroid of unshielded area of vertical 
tail 

&i, k 2 

apparent mass coefficients of body 

L 

lift force 

l 

length of body 


In 

horizontal tail length, distance be- 
tween 0.25 mean aerodynamic chord 
of horizontal tail and airplane center 
of gravity 

Zy 

vertical tail length, distance between 
0.25 mean aerodynamic chord of 
vertical tail and airplane center of 
gravity 

N 

normal force 

NRe 

Reynolds number 

p 

pressure 

Po 

stagnation pressure 

Ps 

static pressure 

Pw 

pressure in wake 

Q 

local dynamic pressure 

Qo 0 

free-stream dynamic pressure, 
2 pUl 

R 

spin radius 

r 

radius 

r c 

corner radius of square cross section 

r 3 

radius of impinging airstream 

r V 

streamwise coordinate of stagnation 
streamline in secondary flow over 
vertical tail 

5 

area 

Sb 

cross-sectional area of body 

"^Bjinax 

maximum cross-sectional area of body 

Ss 

shielded area 

Sv,eff 

effective vertical tail cross-sectional 
area 

U 

local velocity 

U s 

velocity along stagnation streamline 
(fig. 15) 

Uvo 

free-stream velocity 

V B 

volume of body 

X,Y,Z 

axis system 

x,y,z 

coordinate system attached to body 
with origin at center of gravity 

x cg 

center-of-gravity location, measured 
from nose of body 

X C g 

_ £££ 
c 

Xcp 

center of pressure of wing measured 
from leading edge and expressed in 
terms of chord 


3 


w 

2-D 


a angle of attack 

Oil local angle of attack 

/ 3 sideslip angle 

A elemental quantity (strip) 

e correction factor for three-dimensional 

effects over body 

t?h horizontal tail efficiency, qu/Qoo 

Vstall spanwise extent of stall over wing or 

horizontal tail, b ^ r \ 

rjy vertical tail efficiency, ^ 

0 = tan -1 ^ 

% = tan" 1 

A = tan -1 jIj 

= tan" 1 

p density 

0 cross-flow angle, tan” 1 

(j) c cross-flow angle when Cy = 0 

uj reduced spin rate, positive for right 

spin and negative for left spin (as 
viewed from top), Qb/2U 

U angular velocity about spin axis, 

positive for right spin and negative 
for left spin (as viewed from top) 

Subscripts: 

L left 

£ strip or local 

max maximum 

min minimum 

R right 

stall condition or parameter at stall 

Superscript: 

/ parameter in rotational or secondary 

flow effect 

Abbreviations (also used as subscripts): 

B body (fuselage) 

eg center of gravity 

H horizontal tail 

V vertical tail 


wing 

two-dimensional 

Analysis 

Let us consider an airplane model, mounted on 
a rotary balance apparatus in the spin tunnel and 
steadily rotating at a constant angular velocity Q 
about the spin axis, which is vertical. A schematic 
sketch of this configuration is shown in figure 1. For 
the purpose of this analysis, the following assump- 
tions are introduced, mainly because we intend to 
tie the present theory as closely as possible to the 
test conditions employed in the rotary balance tests 
(refs. 11 and 12), which could be generally different 
from those encountered in actual steady-state spin of 
the free airplane: 

1. The spin axis passes through the model cen- 
ter of gravity so that the spin radius R and 
sideslip at center of gravity are zero 

2. The T-axis of the airplane lies in the horizon- 
tal plane so that the bank angle is zero 

3. The effect of power is ignored and the model 
configuration without propeller is considered 
in the analysis 

4. The control surface deflections are zero; in the 
present analysis, the dihedral effects of wing 
and horizontal tail are ignored and the mu- 
tual aerodynamic interferences between vari- 
ous airplane components are ignored except 
the interference between horizontal and verti- 
cal tail surfaces (shielding effect); however, it 
is possible that some of the interference effects 
not considered herein could be significant, as 
discussed later, and need proper modeling 

5. Difference between tunnel test and flight 
Reynolds number is ignored 

In strip theory approach, the surface (say wing) is 
usually divided into a number of strips. The force on 
each strip is assumed to be the same as it would be 
were the surface moving with an equivalent linear 
velocity, that is, the vector sum of all the linear 
velocities at that strip. For a rectangular wing with 
strips oriented normal to the y-axis (fig. 2), this 
equivalent velocity is equal to U + Uy. The strip 
theory approach used for wings and horizontal tails 
is similar to that employed by Glauert (ref. 3). 

Wing Contribution 

Strip theory calculations . Consider the airplane 
model to be rotating in a clockwise (positive) direc- 
tion as viewed from the top. The effect of the angular 
velocity in spin fi is to induce, at any spanwise loca- 
tion y , a velocity component equal to Qy as indicated 


4 



in figure 2. As a result, the local angle of attack a £ 
on the right wing increases and that on the left wing 
decreases compared with the angle of attack at the 
root chord a. The maximum and minimum values of 
a occur, respectively, at the right and left wing tips. 
We have 

ot£ — ot ± 6 (1) 


where 


6 — tan 1 


U 


In equation (1), the + sign applies to the right 
wing and the - sign applies to the left wing; the 
variations in the angle of attack across the wing span 
can be substantial. For the spin-tunnel model of 
reference 14, the figures are as follows: 

For a steep spin, 


ol — 50° <jJ — 0.33 amax — 68.3° a min — 34.3° 


and for a flat spin, 


a = 80° amax - 111.79° a min = 48.21° 


It may be observed that for this case, the flow close 
to the right wing tip is actually reversed; that is, it 
comes from the trailing edge. 

The local dynamic pressure is given by 


Q£ = 




( 2 ) 


The lift and drag forces on any strip of width dy are 

21 


A L = \pUl 


1 + 


and 


A D = \ P Ul 


m- 

(nr 

\U , 


Clx c dy (3) 


Cn £C dy (4) 


Here, Clp and are the sectional lift and drag 
coefficients to be evaluated for the corresponding 
angle of attack ag which depends on the angle of 
attack at the root a and the spin rate Q. 

The normal and chordwise forces on the strip are 
given by 

AN = A L cos(a ± 8) + AD sin(a ± 9) (5) 


(6) 


The aerodynamic coefficients of the wing (based on 
the free-stream dynamic pressure <?oo, wing area S\v, 
and wing mean geometric chord c\y can then be 
derived as follows (refs. 2 and 3): 


C N, W 

l 

2 tan J 

- f [&C N ' R + 

J 0 

sec 4 

e d.e (7) 

C A, W 

1 

2 tan ; 

-Jo [4C ' ui + 

& c A,l\ sec4 

9 de (8) 

C l, W 

1 

4 tan 2 


- 



X tan 9 

sec 4 9 d9 


(9) 

Cn, W 

X 

4 tan 2 

~>Io [AC ^- 

- & c a,fl\ 



X tan 9 

sec 4 9 d9 


(10) 

C m, W 

1 

4 tan 2 

-Jo 

+ A<?AT,l] 



X (x C g - 

- Sep) sec 4 9 dO 


(11) 

where x cp 

is the center-of-pressure location measured 


from the nose of the body (in wing chords). The side- 
force coefficient of the wings Cyw is generally small 
and is neglected: 

Cy, W = 0 ( 12 ) 

In equations (3) to (11), the wing section coeffi- 
cients, C L j and Cj)£i and x cp are input parameters 
which must be known a priori. However, for a partic- 
ular airplane under consideration, if the actual wind- 
tunnel test data are not available, the data given in 
references 19 through 21 can be used. 

Effect of angular velocity on wing . The spin-tunnel 
rotary balance test data (refs. 11 and 12) indicate 
that the aerodynamic coefficients — particularly, nor- 
mal force, rolling, and pitching moments — display a 
pronounced dependence on spin rate. When com- 
pared with these data, the prediction of strip theory 
with static wind-tunnel test data as input differs ap- 
preciably. The reasons for such discrepancies include 
(1) the induced flow from neighboring strips, (2) cross 
flow due to pressure gradients perpendicular to the 
strips, (3) centrifugal forces, and (4) curved flow. 

The effect of spin rate on the aerodynamic forces 
and moments of a spinning wing, as recorded by the 
rotary balance test data (refs. 11 and 12), can be 
summarized as follows: 

Low angles of attack (below stall, say 0° < a < 20°) : 
With an increase in spin rate, 

1. The normal force decreases 


and 

A A = -A L sin(of ± 0) + A D cos(a ± 9) 


5 



2. The rolling moment decreases and becomes 
more damping in nature 

3. The pitching moment increases initially but 
starts dropping subsequently 

High angles of attack (above stall, a > 20° ): With 
an increase in spin rate, 

1. The normal force increases continuously 

2. The rolling moment decreases appreciably and 
becomes more damping; however, for ex- 
tremely high a, the trend reverses and rolling 
moment increases with a 

3. Pitching moment decreases monotonically 

Studies which could provide some physical insight 
into the flow pattern and pressure variation over a 
spinning airplane wing and could have been helpful 
in the development of a mathematical model are al- 
most nonexistent. This is a complex aerodynamic 
problem because it is dominated by extensive regions 
of separated flow, strong viscous-inviscid interaction, 
and configuration and Reynolds number dependence 
in addition to the presence of spin rate. The analysis 
of such a flow problem is considered to be beyond the 
scope of the present study. Instead, we propose to de- 
velop an approximate solution suitable for engineer- 
ing purposes with the help of available experimental 
studies related to this subject. A brief description of 
some related investigations reported in literature is 
presented in the following sections. 

Studies at low angles of attack . In reference 9, 
Himmelskamp has reported wind-tunnel tests on a 
rotating propeller. The model he tested had two 
blades with a chord of 12 cm, a hub diameter of 
10 cm, and a tip diameter of 50 cm (AR = 8.33). (See 
fig. 3.) The cross-sectional shape was profile No. 625 
of the Goettingen collection with a thickness-to- 
chord ratio of 0.20. The advance ratio U/fiR was var- 
ied from 0.124 to 0.459. However, the angle of attack 
of the blade sections was not varied directly; instead 
various combinations of tunnel velocity and advance 
ratio were employed to get the desired angle-of-attack 
variation. Himmelskamp has presented test results in 
the form of sectional lift and drag coefficients based 
on local dynamic pressure at various radial locations 
from 40 to 80 percent radius. These results are given 
in figure 3. The local lift and drag coefficients of the 
blade sections close to the hub are higher than those 
located outboard. Also, the sectional aerodynamic 
lift and drag coefficients are greater than correspond- 
ing two-dimensional static wind-tunnel test data. 
Himmelskamp attributes these increases in sectional 
lift and drag to the phenomenon of local thinning of 
boundary layer near the hub and Coriolis force acting 
favorably on fluid particles. 


If Himmelskamp’s test propeller is considered as 
a rectangular untwisted wing with a chord of 12 cm, 
span of 100 cm, and AR of 8.33 and set in spin- 
ning motion similar to a wing model undergoing ro- 
tary balance tests, we have carried out a strip theory 
analysis using sectional lift and drag data of figure 3 
as input empirical parameters. Calculations extend 
only to u) — 0.35 on account of limitations in data. 
The results of these calculations are presented in fig- 
ure 4 along with the spin-tunnel rotary balance test 
data (ref. 12) for a spinning airplane configuration 
of rectangular wing planform (AR = 6.10). From 
figure 4, there is qualitative agreement despite the 
differences in two geometries. From this, we can in- 
fer the following: 

1. The type of physical flow pattern leading to 
higher values of sectional lift and drag, partic- 
ularly for inboard sections observed on a ro- 
tating propeller (ref. 9) in axial flow, may also 
occur over a spinning wing at low angles of 
attack 

2. With such empirical data which incorporate 
the physics of the flow as input, the strip the- 
ory is capable of providing meaningful esti- 
mates of the aerodynamic coefficients of a ro- 
tating (spinning) wing 

Studies at high angles of attack . The flow 
over a spinning, stalled wing is exceedingly complex 
and is not thoroughly understood. The situation is 
aggravated by a dearth of experimental data which 
treat local flow and pressures. It is necessary, then, 
to deduce what one can from the average or total 
force and moment measurements. 

Application of the strip theory used previously 
for low angles of attack to the stalled wing results 
in the normal- force coefficients presented in figure 5. 
The normal-force coefficient for a = 90° is seen 
to increase by about one third at a; = 1, which 
is consistent with the average increase of dynamic 
pressure due to spinning. Extrapolation of the rotary 
balance test data, however, shows an increase of 
about one in the normal-force coefficient for a = 90° 
and uj = 1. After adjusting for differences at u = 0, 
the discrepancy is seen to grow approximately with 
the square of u). 

Unfortunately, there is no accepted theory which 
accounts for this discrepancy in Cjy. Matters are 
made worse by the manner in which normal force 
is nondimensionalized by dividing by the free-stream 
velocity. A spinning wing at a = 90° will create 
a normal force if the airfoil section is cambered as 
in this paper even when the free-stream velocity ap- 
proaches zero. The consequence is that the nondi- 
mensional normal-force coefficient and cj will grow 
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without bound as U approaches zero. Propeller char- 
acteristics are usually nondimensionalized by using 
the tip velocity which avoids this singularity. 

In reference 13, McCormick hypothesizes that rel- 
atively stagnant air on the lee side of a stalled spin- 
ning wing will have a spanwise pressure gradient due 
to centrifugal force. The resulting reduced pressure 
in the center portion of the wing gives rise to an in- 
crease in the normal-force coefficient (fig. 6); that is, 



where r) s = , the normalized spanwise location of 

°w 

the stalled region. It is assumed that the spin axis is 
at y = 0. The corresponding increment for rolling- 
moment coefficient is 



The expression for C[ w in reference 13 seems to be 
incorrect. 

The corresponding change in the pitching- 
moment coefficient is 

Cm,W = C N , W (®cg “ icp) 

where x(. p is the center of pressure and is assumed to 
be equal to 0.5. 

The addition of the correction term to the normal- 
force coefficient is seen in figure 6 to compare more 
closely with the experimental data. The similar 
comparison in figure 7 for rolling-moment coefficient 
also shows reasonable correlation. Unfortunately, it 
is not possible to validate the detailed assumptions 
integral to the McCormick model because of the lack 
of local flow measurements for spinning stalled wings. 

Body Contribution 

The cross-sectional shape of the body and 
Reynolds number have a significant influence on the 
spin and recovery characteristics of the airplanes 
(refs. 22 through 24). Bodies having rectangular or 
square cross sections with sharp corners particularly 
in the aft portion, produce autorotative side-force 
and yawing moments. These studies indicate that 
rounding off the bottom corners generally produces 
a damping effect whereas rounding off the top comers 
often does not alter the basic prospin tendency. 

Polhamus (ref. 24) has studied the two- 
dimensional axial- and side-force characteristics of 
noncircular cross sections having various corner radii 
at subcritical and supercritical Reynolds numbers. 


He finds that the corner radius and Reynolds number 
have a strong effect on the side-force characteristics 
of two-dimensional noncircular cylinders. In figure 8, 
Polhamus’ experimental data for two-dimensional 
square cylinders of various corner radii are presented. 
For a square section with rounded comers, at sub- 
critical Reynolds numbers, the side force is positive 
(prospin) for small values of <f> and becomes negative 
(antispin) at supercritical Reynolds numbers. How- 
ever, for the section having an 8-percent corner ra- 
dius, the Reynolds number has very little effect on 
the side-force variation. This pronounced effect of 
the Reynolds number may sometimes cause the spin- 
tunnel free-flight model to exhibit a flat spin mode, 
whereas the full-scale airplane may refuse to enter a 
flat spin (ref. 15). Such complex behavior of noncir- 
cular cross-sectional bodies typical of airplane bodies 
requires considerable care in making theoretical anal- 
yses of spinning airplanes. 

During a spin, the body operates at a high angle 
of attack, and its cross sections experience a vary- 
ing sideslip. At present, no comprehensive method 
is available for the prediction of aerodynamic char- 
acteristics of noncircular cross-sectional bodies oper- 
ating at combined high angles of attack and sideslip. 
Therefore, in the following sections, a semiempirical 
procedure is developed for a static body operating 
at (1) angle of attack, (2) sideslip, and (3) combined 
angles of attack and sideslip. Then this procedure is 
extended to the problem of a spinning body with the 
use of strip theory. 

For the analysis of static body at angles of attack 
and sideslip, a coordinate system as shown in figure 9 
is assumed where the origin is the location at the 
nose. However, for the spinning body, the same 
coordinate system as defined earlier with the origin 
at the center of gravity (fig. 2) is employed. 


Body at angle of attack. According to Allen 
(ref. 25), for a body of revolution at small angle of 
attack, 


Cl, b = 


(k 2 — ki ) sin 2a cos 


a. rl 


Sb,v 


fdS b 

Jo 


dx 


+ 2 sin 2 a cos a f l 
%,ma x Jq 


ercd dx 


(19) 


Cd, b = 


(&2 — ) sin 2a sin j f l dS& 


+ 


2 sin 3 a f l 

t,max Jq 


f ds b 

Jo dx 


dx 


Sb 


ercd dx 


( 20 ) 
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follows: 


(fc 2 -fci)sin2acosf f l dS B , 

C m ,° jr 1 1 -2- <*, - ») * 

‘ 2c '^ e a J (xcg - x) dx (21) 

In equations (19) to (21), the first term represents the 
potential term proposed originally by Munk (ref. 26), 
and the second term is Allen’s modification (ref. 25) 
to Munk’s theory derived from the viscous cross-flow 
concepts containing the empirical two-dimensional 
drag coefficient (e^) of the circular cylinder. The 
zero-lift drag term appearing in Allen’s equation 
(ref. 25) is ignored. Jorgensen (ref. 27) extends this 
approach to noncircular bodies with the following 
modification: 


C L, b = 


C\ (k 2 — ^i) sin 2a cos 


Sb,e 


jr f‘ dS s 

Jo dx 


dx 


sin 2 a cos a f / 00 \ 

-5 / tb h C x , c d* (22) 

°B,max J o 


C D, B = 


Ci (At 2 — &i) sin 2a sin 


S B,n 


’ n § f l d&B 

Jo dx 


dx 


[ l (bh C XtC dx (23) 

^E^max Jo 


Cl (fc 2 -All) sin 2a cos § f l dS B , 

C m,B = J Q (*cg - x) dx 

+ J tb h C XiC (x C g - x) dx (24) 

The modifications are (1) the first term is mul- 
tiplied by a constant C\ so that the quantity 
Ci (&2 — ki) represents the apparent mass coefficient 
of the noncircular cross section and (2) the coeffi- 
cient Cx, c (based on 6/J replaces c j and represents 
the cross- flow drag coefficient of the noncircular cross 
section. It may be noted that for a circular cross sec- 
tion, Ci = 1.0, bh — 2 r, and Cx, c = c d- 

Jorgensen (ref. 27) proposes that equations (22) 
to (24) be applied for a = 0° to 90°. In the present 
analysis, all the concepts proposed by Jorgensen 
(ref. 27) are used. 

If Cl and Cp are known, the normal- and 
axial-force coefficients can be determined from equa- 
tions (5) and (6). Since the axial force does not 
produce any moment about the spin axis, it is not 
addressed. 


Cy,B = 


Cl (k 2 -k 1 )sin 2/? cos Q f l 




f 

JO 


dS'Q 

dx 


+ 


sin 2 /? 

■^B, max 



dx 


(25) 


and 

Ci(k 2 -fci)sin2/?cosf f l dS B , 

c„, B = ^ J 0 (Icg ~ x) dx 

• 2 a fl 

+ J h (C Y,c (leg - x) dx (26) 

In equations (25) and (26), notice that the first term 
has a negative sign and the empirical side-force coef- 
ficient Cy c (which is supposed to embody the side- 
force characteristics of a noncircular cross section) 
replaces Cx, c ■ In equations (22) through (26), C Xc 
and Cy c are empirical parameters. Therefore, to 
proceed with the calculations, it is necessary to have 
this information from experimental data. 

The rolling moment developed by the body is 
generally small because of its limited moment arm 
and is ignored. 


Combined angles of attack and sideslip. For any 

combination of a and 3 (o f 0), we can define a 
cross-flow angle <f> (fig. 9) such that, 


tan <j> — 


tan/? 
sin a 


(27) 


Equation (27) relates the sectional cross-flow angle (j) 
with angles of attack and sideslip. To evaluate the 
aerodynamic forces and moments under combined a 
and /?, the proposed procedure is as follows: 

1. For the given values of a and /?, find the value 
of <f> from equation (27) 

2. Obtain the corresponding values of Cx, c and 
Cy, c f° r the subject cross-sectional shape from 
the experimental data which, as said before, is 
supposed to be available 

3. With these values of Cx iC and Cy ?c , solve 
equations (19) through (27) to obtain the re- 
quired aerodynamic coefficients 


Strip theory calculations . The variation of cross- 
flow angle and dynamic pressure along the length of 
a rotating body are given by 


Body in sideslip . For a noncircular body in 

sideslip, we extend Jorgensen’s concept (ref. 27) as 


<t>(x) = tan 1 


Qx 

Uoo 


(28) 
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(29) 

Notice that the cross-flow angle <f> is independent of 
the angle of attack. 

The strip theory procedure for the evaluation of 
aerodynamic characteristics of a rotating (spinning) 
body is as follows: 

1. Divide the subject body into a convenient 
number of axial strips (say 25) of length dx 

2. Then at a given strip, evaluate <j>{x) and 
Qi,b( x ) f rom equations (28) and (29) 

3. From the following equations, evaluate the 
aerodynamic coefficients (per emit length) of 
each strip: 


ACw b = 


Ci (&2 - fci) sin 2a cos % dS& 


sin 2 a , 

+ c eb h.Cx,c 

^B,max 


ACvb = 


— AC^B y ^ ^Cg X) 


(31) 

C\ (&2 - k \ ) sin 2/d cos | 

r dSft 


^B,max 

dx 


+ f S <hCy,. 


(32) 

a s>t *^B,max / \ 

= Acy >B y (icg x) 


(33) 


4. Numerically integrate the following expres- 
sions to obtain the aerodynamic coefficients of 
a rotating (spinning) body: 


f l 1 I", /flzsinaX 2 

Cn - b= U . 

x A Cn,b( x ) dx 
_ f ll [, /flxsina\ 2 

Cm ’ B ~J- h I; + \~u^T ) . 

X AC m , B (x) dx 
_ , f l i |\ /flxsina\ 2 

c ™=L [ i+ {— )J 

X AC r y ? B(x) dx 

_ f l * L /fix sin a\ 2 


x AC„ iB (i) dx (37) 

C/,B = 0 (38) 

where l\ and 1 2 are the distances of the center 
of gravity from leading and trailing edges, 
respectively (fig. 10). 

Rotational flow effects. Clarkson, Malcolm, and 
Chapman (ref. 28) have investigated the effect of an- 
gular velocity fl on the aerodynamic forces and mo- 
ments of certain autorotative flow mechanisms, one 
of which is the body with rectangular cross sections 
investigated by Polhamus (ref. 24). In their tests, the 
Reynolds numbers ranged from subcritical to super- 
critical, and the nondimensional spin rates based on 
the body length varied from —0.2 to 0.20. They ob- 
served that the autorotative flow mechanisms, once 
formed, get locked in and rotate with the body. The 
speed of rotation fl apparently did not have a signif- 
icant influence on the force and moment characteris- 
tics. Although the range of fl considered in this work 
is much higher, for the lack of more comprehensive 
information, the rotational flow effects are ignored 
for the body. 

Horizontal Tail Contribution 

Strip theory calculations. Proceeding similar to 
the analysis of the wing, we obtain 


c "' h = L " {acn ' kr + ac ' v ' hi) 

x sec 4 #h dOft 


Ca » - 2Sw”S ^ /„ " (AC ' 1 ' HR + AC * Ht) 

x sec 4 Oft ddft (40) 

Cm,H = ~ C N, H (W g w) ( 41 ) 


x tan 0 h sec 4 % dOn 

r = m gH&H [ Xh 
n,H 4 tan 2 Ah Syi/byf Jo 

x tan On sec 4 Oft dOft 


dOft (42) 

[ (&Ca ,hr ~ ACa,hl) 

Jo 

16ft (43) 


In equations (39) through (43), ?7 h is the horizon- 
tal tail efficiency. We assume that rjft = 1 if the hori- 
zontal tail is clearly out of wing wake and rfn = 0 if it 
is fully immersed in it (fig. 10). The wake of a stalled 
wing can be assumed to be bounded approximately 
by the straight lines extending from the leading and 
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trailing edges of the wing in the direction of the local 
flow (ref. 19). 

Rotational flow effects. For the rotational flow 
effects on the horizontal tail, the procedure is similar 
to th at for the wing calculations except the parame- 

ter h replaces y for the distance of the strip 

from the center of rotation. In this case, the spanwise 
pressure variations leading to the rotational flow cor- 
rections are the ^-components of the radial pressure 
gradient term. With this, we obtain the following re- 
lations for the horizontal tail with the z-component 
pressure forces arising because of the assumed radial 
pressure gradient being ignored: 

(44) 

= 0 ( 45 ) 

and for partial separation on either panel, 

2 c* 

C N,U ~ H (l + ^stall) ( 46 ) 

W _ W 2 7?h-5h% (■, „4 

- ~ 165 w i-w V “ (47) 

C'm,H= -C'n,h(^) (48) 

Vertical Tail Contribution 

Shielding effect. When stalled, the flow over the 
horizontal tail breaks off from its upper surface. This 
separated flow forms a region of dead air or wake on 
the lee side of the horizontal tail. On a spinning 
airplane, the wake flow of the outboard panel (the 
left panel in a right spin and the right panel in a 
left spin) of the horizontal tail is pushed towards the 
vertical tail (fig. 11). This results in the shielding of 
the vertical tail by the stalled stagnant fluid of the 
wake. Notice that the wake of the inboard panel is 
swept away and does not create any shielding over 
the vertical tail surface. Thus, only the outboard 
panel of the horizontal tail would be responsible for 
the shielding phenomenon. 

In order to determine the aerodynamic forces and 
moments over the shielded vertical tail, the following 
assumptions are introduced: 

1. The shielded area is a function of angle of 
attack and is independent of angular velocity 
in spin 


2. The characteristics of the wake of the horizon- 
tal stabilizer can be approximated on the basis 
of experimental data on the two-dimensional 
flat-plate wake 

3. The unshielded region of the vertical tail re- 
tains its effectiveness 

4. The shielded area of the vertical tail is other- 
wise ineffective aerodynamically, except for 
the secondary flow effects discussed later 

Modeling of two-dimensional flat-plate wake. At 

low Reynolds numbers, the mean flow in the wake, 
close to the flat plate, may be broadly divided into 
two regions (fig. 12(a)): (1) a recirculating region 
(ARB) bounded by the flat plate and the dividing 
streamline and (2) the external, irrotational flow. 
(See ref. 29.) 

Arie and Rouse (ref. 30) find that the recircu- 
lation region (ARB) has an approximately ellipti- 
cal shape and extends to a downstream distance of 
about 2.4 times the flat-plate chord (cy). Abernathy 
(ref. 31) gives the following relations for a\ and df 
(fig. 12(a)) of the recirculating region, 

ai = 0.5cy 

df = \pic j sin a (49) 

If these two sources of experimental data (refs. 30 
and 31) are combined, the part of the wake behind 
the normal flat plate can be sketched as shown in 
figure 12(a). The major axis a of the ellipse turns out 
to be equal to 3.8 cy, and the part of the ellipse ahead 
of the flat plate in the upstream direction (which is 
of no physical significance) has a length of 1.4cy. For 
application of this result to other angles of attack, 
we assume that the correlation with angle of attack 
given in equation (49) holds for other dimensions of 
the ellipse also. In other words, 

a = 3.8cysina (50) 

ai = 0.5cy sin a (51) 

Thus for any given a, the geometrical parameters 
such as a, a i, and dr can be determined from equa- 
tions (49) through (51), and the shape of the ellip- 
tical recirculating bubble can be sketched as shown 
schematically in figure 12(b). The shielded area ( S s ) 
of the vertical tail corresponds to the shaded region. 

The experimental results of Abernathy (ref. 31) 
based on which model is constructed are obtained 
for a > 40°. Therefore, for a sta u < a < 40°, we 
use linear interpolation with the assumption that at 
a = c*stalb we get S s = 0. 
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Calculation of aerodynamic coefficients. Here, 
a direct approach can be made instead of the strip 
theory because the vertical tail normally has a small 
aspect ratio. The sideslip and dynamic pressure at 
the vertical tail are given by 


0y = - sin 1 < 

2u; (/y/6\y) sin a 1 

(52) 

\/l + [2a; flv/6\v) s * n ) 

l 


1 ,,2 

QV = ^P U 5oW 

1 ! / 2 ujly sin a \ 2 

V 6 W / 

(53) 


Under a partially shielded condition, we assume that 
the side force is given by 


Fyy = QV C Y,V S V,efi 


f Y,V = \pUloVV 




(54) 

Cy.V^V.eff (55) 


where 5y >e fr = Sy — S s and Cyy is the side-force 
coefficient of the vertical tail. Its determination is 
discussed later. Then 


C Y,V 


f y,v 
q s w 


(56) 



c iy 


n h v 
c r,vr- 
6 W 


(58) 


C n,V 


r v 

-Cy,vr 


(59) 


where hy is the vertical distance between the 
X-axis and the centroid of the unshielded area (taken 
positive upwards) of the vertical tail. In left spin, to 
is negative, 0y is positive, and Cy, v is negative. As 
a result, C/ y and C n y have opposite signs to those 
given in equations (58) and (59). Further, we assume 


CjV.V = C m,V = Ca,V = 0 (60) 


Secondary flow effects . The spin-tunnel rotary 
balance data (ref. 11) indicate the existence of a no- 
ticeable amount of damping in yaw (in excess of that 
predicted by eqs. (52) through (54)) at combined 
high angles of attack and spin rates. Such a fact, 
however, may never be noticed by the pilot in spin 
flight tests because the magnitude of such a damp- 
ing effect may be small under spin flight test condi- 
tions. This development of the damping in yaw is an 
interesting phenomenon because it occurs when the 


vertical tail is assumed to have lost its effectiveness 
due to shielding. A detailed analysis of this prob- 
lem which is characterized by extensive regions of 
separated and recirculating flows, in addition to the 
presence of the spin rate, is considered to be beyond 
the scope of this study. However, a physical expla- 
nation and an approximate mathematical model of 
this phenomenon, referred to here as the “secondary 
effect,” is presented in the following discussion. 

A simple way of visualizing this secondary effect is 
to assume that the airplane model is in a pure rotary 
motion, and the upward velocity in the spin tunnel 
is zero. Here, the “flat-plate” drag of the vertical 
tail generates a positive rolling moment and a nega- 
tive (damping) yawing moment about the center-of- 
gravity of the airplane model. This phenomenon is 
analogous to the one referred to as secondary effect. 
The difference here is that the static pressure in the 
air surrounding the vertical tail is equal to p w , the 
pressure in the wake of the horizontal tail. 

On a spinning airplane, the mass of fluid in 
the wake has a relative transverse velocity equal to 
H/y sin a with respect to the airplane and is assumed 
to impinge perpendicularly on the windward side of 
the vertical tail. As a result, a stagnation line forms 
over the windward side of the vertical tail about 
which the airstream divides itself and moves in oppo- 
site directions (fig. 13). The existence of such a flow 
pattern has been noticed in the full-scale flight tests 
conducted by the Langley Research Center. 

As a consequence of this wake flow impingement, 
which is designated here as “secondary flow effect,” 
a side force Fyy (positive in right spin and nega- 
tive in left spin) develops on the fin and generates an 
autorotative rolling moment and an antispin yawing 
moment (fig. 14). It is assumed that this wake flow 
impingement is limited to the shielded portion of the 
vertical tail. An estimation of the incremental aero- 
dynamic forces and moments due to the secondary 
flow effect on the vertical tail is done as follows. 

It is assumed that the wake flow impingement 
pattern is similar to that of the impingement of a 
uniform flow of radius rj (fig. 15). (See ref. 32). For 
uniform flow impingement on a normal surface, the 
pressure and velocity along the stagnation streamline 
are related by Bernouli’s equation, 

Pa + ^ pUs = Po (61) 

where p 9 ,U s , and p 0 are static pressure, velocity, and 
stagnation pressure. Also, 


Po = P + ^ pU 2 


(62) 
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where p and U are the pressure and velocity of the 
uniform impingement stream. Close to the stagna- 
tion point, we can assume (ref. 32) 

Us = -U (63) 

r 3 


For a partially shielded condition, replace Sy by S s 
to give 


w _ USs fwly sin a \ 2 

°yy - 35w {~~hT) (71) 


where r is the distance from the stagnation point 
measured along the surface streamline, and ry is 
the radius of the impinging airstream. Substituting 
equations (62) and (63) into equation (61) results 
in the following equation, which gives the pressure 
distribution on the surface, 


Ps = P + ^ pU 2 



(64) 


Now for application of this model to wake flow im- 
pingement over the vertical tail, we have (fig. 13) 


U = Q/y sin a ' 
P = Pw 
r j = C V 
r = ry 


(65) 


Therefore, the pressure distribution along the chord 
of the vertical tail is given as 


Ps = Pw + 2 P (Mv sin a ) 2 



( 66 ) 


On the lee side of the vertical tail, pressure is as- 
sumed to be equal to p w . Thus, the net side force on 
a strip with width Ary due to the pressure difference 
(Pa ~Pw) is given by 

A Fy V = (p. - p w ) by Ary (67) 

A Fyy = i p(CUysina ) 2 [l - ( ry/cy ) 2 ] by Ary (68) 

The total side-force coefficient on the vertical tail 
when it is fully shielded is 


Then the rolling- and yawing-moment coefficients 
due to this secondary flow effect are (fig. 9(b)) 


C 'n,V = TC Y,V]^ J 


(72) 


For right spins, C\ y is positive and C? n y is negative; 
for left spins, C[ v is negative and C f n y is positive. 
Thus, from equations (70) and (72), we observe that 
the magnitudes of the side-force, rolling-moment, 
and yawing-moment coefficients due to the secondary 
flow effect on the vertical tail will assume significance 
at combined large angles of attack and spin rates. 


Total Airplane 

Expressions (eqs. (73) through (80)) are presented 
for determining the aerodynamic coefficients of var- 
ious combinations of airplane components based on 
the calculations in the previous section. The format 
for the presentation of results follows the manner in 
which the rotary balance test data are reported in 
references 11 and 12. It may be noted that the axial- 
force data are not available in reference 11. Since the 
axial force does not play any significant role in equi- 
librium spin, it is not considered. However, for com- 
plete airplane configurations, the axial-force data are 
available in reference 12. Therefore, in the present 
paper, results on the axial force are limited to air- 
plane configurations only. The prime (/) denotes the 
rotational or secondary flow effect. 


_ 2by {Qly $ina\ 2 f°y/ 2 

Cy * ~ % u~) Jo 



lf _ llSy /u;/ysina\ 

Y,v 3Sw V &w / 


Body B configuration: 



Cn = CV, B 

(69) 

03 

£ 

II 

£ 


C m — Cj 

(70) 

Cn = C n> B 


( 73 ) 
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BW configuration: 


BWV configuration: 


C N = B + Cjv, w + C' N w ' 
Cy = Cy, B 

Cm ~ C m> B + C m ,W + Cm,W f 

Q = Q,w + C?,w 

Cn = Cn,B + Cn^w 


(74) 


Cy - Cy jB + CjV.W + Cy W 
Cy = Cy )B + Cy,y 
C m — Cm,B + ^m,W + C’m.W ' 
Cj = Q,w + c?,w + Q,v 
Cn = C nj B + C Wj w + C n y , 
BHV configuration: 


( 78 ) 


BH configuration: 


C N = Cy )B + C^h + Cy jH 
Cy = Cy >B 

Cm — C m B + C m ,H + C^ H ► 
C l = Q,H + Cj H 

Cn - C n B + C n>H 


(75) 


BY configuration: 


Cn = Cy >B 
Cy = Cy B + Cyy 
Cm = C mB 

Cl = Ciy 

Cn — C n? B 4" C n? y 


(76) 


BWH configuration: 


Cat = Cy,B + Cy ? w 4 Cy,w 4- Cy,H 4- C^ H 
Cy = Cye 


Cm — C m B 4 C m> w + C^ vV + C m? H + C^ h > 

C; = Q,W 4- C/ ? w + C/ } h + C/ H I 


(77) 


Cn = C n? B + C n> w 4 C n> H 


Cat = Cy B + Cy 5 H + Cy H 
Cy = Cy ? B + Cyy + Cyy 
Cm = C m 3 + C m? H + C^ H 
Cl — Q,H + C/ H + C/ ? y + Cj y 

Cn = C n)B + C n?H + C nj y + C^y 


(79) 


BWHV configuration: 


Cy = Cy ?B + Cy^w + Cy\y 4 Cy iB + Cy^H 
c^ = C>i ? b + Ca,w + Ca,h 

Cy — Cy )B 4 Cyw + Cyj| 4 Cyy 4 Cyy 
Cm = C m B 4 C m , W + C^y/ + C m ,H 4 - Cm,H 
Ci = C^w 4 C/ w + C/ 5 H + C/ H 4 C/y 4 - c/y 
C n = C n?B 4 C n> w + C n> H + C n y 4 C^ y 


(80) 


Application 

The semiempirical aerodynamic analysis of the 
spinning airplane developed herein has been applied 
to a light, general aviation, spin research airplane, 
called “model A,” which is extensively studied in the 
stall/spin research program of the Langley Research 
Center (refs. 9 through 13). 

The dimensional characteristics of the full-scale 
model A airplane are presented in table I. A three- 
view drawing of the basic airplane is shown in fig- 
ure 16, and the details of tails 3, 4, and 5 for which 
the calculations are performed are given in figure 17. 
The tail configuration of the basic airplane is des- 
ignated as tail 4, and wing section is the modified 
NACA 642-415 airfoil. The major difference in vari- 
ous tail configurations is the location of the horizon- 
tal tail with respect to the vertical tail. Among these 
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configurations, tail 3 produces a maximum shielding 
effect, and tail 5 (T-tail) is not expected to produce 
any shielding. These tail configurations are selected 
to examine the validity of the shielding model devel- 
oped earlier. 

Input Data 

In the semiempirical theory of a spinning airplane 
presented, several empirical constants are embedded. 
In the following discussion, these constants are eval- 
uated for the subject airplane model. 

Wing. The empirical constants are Clj, Cpj, 
(or Cpjj and C,y) and z cp - In reference 10, the 
static wind-tunnel test data are presented for the 
body B and the wing-body BW combinations. Ig- 
noring wing-body interference, the wing characteris- 
tics are derived by subtracting the data for the body 
from the data for the wing-body combination. The 
experimental center-of-pressure z cp data are taken 
from reference 33. These derived data are presented 
in figure 18. 

Body. The empirical constants to be evalu- 
ated for the subject body are Ci, Cx, C i and Cy, c - 
The constant C\ depends on the body length and 
cross-sectional shape, and Cx jC and Cy, c are aero- 
dynamic parameters which depend on the cross- 
sectional shape, cross-flow angle <£, and the Reynolds 
number. The range of cross-flow angle <j> over which 
the experimental data for Cx tC and Cy c are required 
is shown in figure 19. For rotary (spin) conditions, 
the maximum value of <t> is around 40° for all angles 
of attack. However, in static calculations, larger val- 
ues of <f> (up to approximately 80°) are encountered. 

The coefficients Cx, c and CV,c correspond to a 
two-dimensional cylinder which has a constant cross 
section. However, the actual cross-sectional shape of 
the subject body varies along the length as shown in 
figure 20. Therefore, we have to idealize the body 
of model A, as the one having the same values of 

and dS^/dx as the subject body (fig. 21) but a 
constant cross-sectional shape. For this purpose, we 
refer to the study of Bihrle and Bowman (ref. 23) who 
tested various bodies similar to that of the model A 
airplane. In their tests, the aft body shape was fixed, 
and the cross-sectional shape of the forward part 
was varied. They observed that the cross-sectional 
shape of the nose did not have significant influence on 
the autorotational characteristics of the body. This 
observation may be attributed to the short moment 
arm of the forward part of the body. However, these 
comments are specific to this configuration and may 
not be valid for other bodies. Thus, we assume that 


the idealized constant cross-sectional shape of the 
subject body is that of the aft portion, which happens 
to be a square section with sharp bottom corners 
and a well-rounded top surface as shown in figure 22. 
For this idealized section, the empirical constants are 
evaluated as follows: 

C\: For a square section with sharp corners, 
Jorgensen (ref. 27) suggests a value of 1.19, 
which has been used in the present calcula- 
tions. 

Cx,c : Wind-tunnel test data are not available for the 
idealized cross-sectional shape. Perhaps the 
closest data available are those reported by 
Polhamus (ref. 24) for a square section with a 
corner radius of 8 percent (fig. 8). With the as- 
sumption that the nature of windward corners 
has the primary influence on the cross-flow 
characteristics compared with lee-side comers, 
we use these data in the present calculations. 

Cy c : The required test data are not available. The 
generation of these data is discussed in the 
following paragraphs: 

Side force — In reference 10, the three- 
dimensional static side-force data on the iso- 
lated body of the model A airplane are pre- 
sented at various angles of attack between 0° 
to 90° and sideslip angle of 0° to 30°. With the 
use of equation (27), the variation of side force 
with cross-flow angle (j> is obtained as shown in 
figure 23. The test data corresponding to high 
angle of attack (a > 70°) are well correlated 
because the cross flow is very much indepen- 
dent of the axial flow at these angles of attack. 
However, the large scatter in the data points 
corresponding to tests at lower angles of at- 
tack is because of the strong interdependence 
between axial and cross flows. Therefore, a 
meanline is drawn through the scattered data 
to approximate the variation of Cy jC at high a 
with (f> as shown in figure 23. These correlated 
test data are designated as “Cy )C (B).” 

Square cylinder test data — A cross-sectional 
shape which is close to that of an idealized 
body section and for which comprehensive 
two-dimensional test data are available is the 
square section having an 8-percent corner ra- 
dius (ref. 24), which has been discussed ear- 
lier. The values of Cx, c and Cy c for this body 
are nearly independent of the Reynolds num- 
ber (ref. 24). However, for the present ideal- 
ized cross section, the cross-flow aerodynamic 
forces can have a dependence on the Reynolds 
number on account of the rounded top surface 
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(fig. 24). These experimental data are desig- 
nated as “Cy, c (SC).” 

These two sets of data are used in the present 
calculations as follows: 

Static aerodynamic coefficients: 

a < 70°, C Xc and Cy c (SC) 

> 70°, C x \ c and C Y ] c (B) 

Rotary aerodynamics: 

All angles of attack, C x , c an d Cy ?c ( B) 

Horizontal tail . The empirical constants required 
are C/v,h, an d z C p,H* These coefficients are 

assumed to be identical to those of the wing as given 
in figure 18. It was assumed that the differences in 
the aspect ratio and airfoil shape do not significantly 
affect these parameters. The determination of 77 h, 
the tail efficiency parameter, is determined graphi- 
cally by considering whether the tail is submerged 
(r?H = 0) or falls out of the wing wake (??h = 1) as 
illustrated in figure 25. Tails 3 and 4 are out of the 
wing wake for a > 20°. However, the T-tail (tail 5) 
gets immersed in the wing wake around a = 20° and 
can be assumed to fall out of wake for a > 30°. 

Vertical tail . The empirical constant Cy\ ap- 
pearing in equations (54), (55), and (57) can be eval- 
uated with Datcom methods (ref. 19) for sideslip an- 
gles below stall. Beyond stall, assuming that the low 
aspect ratio vertical tail behaves like a square plate, 
Hoerner’s data (ref. 20) can be used for the verti- 
cal tail operating beyond stall in sideslip. Based on 
these two sources of data, the side-force data for the 
vertical tail are generated as shown in figure 26. 

Figure 27 presents the schematics of the graphical 
determination of the shielded area of the vertical tail 
at various angles of attack. The calculated variation 
of S s with a for tails 3 and 4 is shown in figure 28. 
It may be noted that S s = 0 for tail 5 (T-tail 
configuration). 

Presentation of Results 

The static longitudinal coefficients of the body 
(B) are shown in figure 29, and the lateral aerody- 
namic coefficients are shown in figures 30 and 31. 
The static side-force and yawing-moment character- 
istics of the wing-body (BW) combination at selected 
angles of attack are given in figure 32. The predicted 
static side-force and yawing-moment coefficients of 
V and BHV configurations, which are based on the 
present shielding model, are given in figure 33 along 
with static wind-tunnel data (ref. 10). 


The steady-state rotary (spin) aerodynamic char- 
acteristics of basic airplane components are presented 
in tables II through VI. With these values and equa- 
tions (73) through (80), the aerodynamic coefficients 
for any desired configuration can be obtained. These 
results are plotted in figures 34 to 44, along with 
corresponding spin-tunnel rotary balance test data 
(refs. 11 and 12). Such a detailed, piece-by-piece 
comparison helps to bring out the capabilities and 
limitations of the present approach. In figures 45 
to 49, the aerodynamic coefficients are presented 
for different tail configurations and at angles of at- 
tack around which equilibrium spin modes have been 
recorded in free spinning tests in the spin tunnel 
(ref. 14). These results are also compared with 
the corresponding spin-tunnel rotary balance data 
(ref. 12). 

Discussion 

Static Aerodynamic Characteristics 

This discussion is limited to body and vertical tail 
because the static aerodynamic coefficients for other 
components (wing and horizontal tail) are directly 
taken from the wind-tunnel test data (ref. 10) and 
not calculated here. 

Body. The predicted value of Cd (fig. 29) agrees 
well with static wind-tunnel test data up to a = 60°, 
but for a > 60°, the present value is on the higher 
side. At a = 90°, the difference between the two 
results is about 15 percent. 

C^. Predicted values of Cl are in fair agreement 
with experimental data (fig. 29) for low (a < 20°) 
and high (a > 70°) angles of attack. For 20° < a < 
70°, the significant differences exist. The maximum 
value of Cl and the angle of attack at which this 
occurs are fairly well predicted. 

C m * The basic static stability trend is captured 
by the present approach; that is, the static instability 
up to a = 30° and subsequent stability are predicted. 
However, significant differences in magnitudes are 
noted to exist. 

Cy. The side-force coefficient Cy has an inter- 
esting variation with a and 0. At low angles of at- 
tack, for example, a < 20°, the side force is negative 
at all values of 0 (fig. 30) and is reasonably well- 
predicted by the present theory. At higher angles of 
attack with the exceptions of 30° and 40°, the side 
force is positive and changes sign for higher values 
of 0. Although calculations pick up this trend well, 
significant differences exist between theory and ex- 
perimental data. 
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C n . At all angles of attack, C n is generally nega- 
tive for the range of /? shown in figure 31. According 
to static wind-tunnel data, the body is directionally 
unstable up to /? = 15° to 20° and becomes stable 
subsequently for higher values of /?. This trend is 
captured by the present method, particularly at high 
angles of attack. As before, significant differences in 
magnitudes exist. 

BW configuration . In figure 32, the present cal- 
culations (ignoring wing-body interference effect) are 
compared with static wind-tunnel test data (ref. 10). 
Also included in this figure are the experimental data 
for the body (B) for relative comparison. It is ob- 
served that the significant interference effect exists 
between wing and body. 

BHV configuration . From figures 33 and 34, 
we observe that the predicted values of Cy based 
on the shielding effect model are in fair agreement 
with static wind-tunnel test data (ref. 10). However, 
the calculated values of C n differ considerably from 
the static wind-tunnel data (ref. 10) particularly at 
high angles of attack. At these angles of attack, 
the vertical tail has lost most of its effectiveness. 
Therefore, this discrepancy in C n could be due to 
the aft body-tail interference effect (not modeled 
here) which may alter the side force locally without 
affecting overall side force much but changes the 
yawing moment significantly on account of its large 
moment arm. 

Rotary Aerodynamic Characteristics 

Body (B) configuration . From figure 34, it is 
observed that large differences exist in magnitudes 
between the present result and rotary balance test 
data (ref. 11). Although the basic autorot at ional 
tendency ( C n > 0 for u > 0) of the body is predicted, 
the autorotational speeds (value of uj when C n — 0) 
are not well predicted. 

BW configuration . The approach of strip theory 
combined with McCormick’s rotational flow model 
(ref. 13) gives a good result in comparison with rotary 
balance test data for <7#, Q, and C m (fig. 35). The 
increments to Cjy, Q, and C m calculated on the 
basis of McCormick’s rotational flow model are of 
considerable magnitude at high spin rates. However, 
significant differences between the present results and 
rotary balance test data exist, particularly for side- 
force and yawing-moment coefficients. 

BH configuration . The variation of Cy , Cjy, 
and C m follows the earlier trends (fig. 36). For 


the BH configuration, the calculated rolling-moment 
coefficient is very small (e.g., Q « 0 for a = 
90°). However, the rotary balance data indicate the 
existence of substantial rolling moment (autorotative 
type) at a = 60° and 90°, which is far in excess of 
that given by the strip theory. 

In figure 36, for the coefficients Cy and C n , the 
rotary balance test data for the body alone are also 
superposed. Observe that in relative comparison, the 
yawing-moment coefficient for the BH configuration 
is significantly different from the body-alone values. 
The strip theory contribution of the horizontal tail to 
Cy and C n is very small. Therefore, the significant 
change in C n when both the body and horizontal 
tail are together must be on account of the aft body- 
horizontal tail interference as noted by earlier authors 
(refs. 22 and 23) and not considered here. 

BV configuration . In figure 37, the discrepancies 
in Ci and C n for the BV configuration are less severe 
compared with those for the BH configuration. The 
interference effect between the aft body-vertical tail 
is also significant as seen by a relative comparison for 
BV and B data for Cy and C n . Although the main 
trend of results is consistent with rotary balance 
tests, differences in terms of magnitude still exist. 

BWH configuration . The coefficients Cj y, C/, and 
C m (fig- 38) follow the trends noted earlier. However, 
the predicted side-force and yawing-moment coeffi- 
cients differ appreciably from the rotary balance test 
data. Contributions of wing and horizontal tails to 
side-force and yawing-moment coefficients are small 
in comparison to that of the body. 

In figure 38, the rotary balance test data for 
Cy and C n for B, BW, and BH configurations are 
also included along with the data for the BWH 
configuration. From a relative comparison of these 
values, it is clear that significant interference exists 
between wing, body, and horizontal tail surfaces. 
As stated earlier, these interference effects are not 
modeled here. The magnitudes of interference effects 
are so appreciable that the presence of wing and 
horizontal tails has masked the basic autorotational 
tendency of the body. 

BWV configuration . The basic trends of most of 
the experimental rotary aerodynamic coefficients are 
predicted by the present method (fig. 39). Differences 
in magnitudes do exist. The disagreement in Cy 
suggests that there is an appreciable interference 
effect (between body, wing, and vertical tail) as noted 
previously. 
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BHV configuration . The predicted trends of Cjy, 
C m , and C[ agree with those given by rotary balance 
data (figs. 40 and 41). The calculated side-force 
and yawing-moment coefficients differ considerably 
from the rotary balance test data. This phenomenon 
is attributed to the aft body-tail interference effect 
noted earlier. 

BWHV configuration . The present theory pre- 
dicts the characteristics of the total airplane reason- 
ably well at low and moderate angles of attack and 
low spin rates. The differences in magnitudes be- 
come significant at high angles of attack and spin 
rates, particularly for side-force, rolling-, and yawing- 
moment coefficients. 

Prediction of Equilibrium Spin Modes 

For the spin research airplane model A, the free 
spinning model tests (ref. 14) have indicated the 
existence of steady-state (equilibrium) spin modes as 
follows: 

Tail 3: 

(1) Moderately steep spin mode, a = 50° and 
u) = 0.33 

(2) Flat spin mode, a = 80° and u = 0.62 
Tail 4: 

(1) Moderately steep spin mode, a = 35° and 
w = 0.22 

(2) Flat spin mode, a = 77° and oj = 0.92 
Tail 5: 

Moderately steep spin mode, a = 41° and 

u = 0.234 

In figures 45 to 49, the predicted rotary aerody- 
namic coefficients for these spin modes are presented 
and are compared with corresponding spin-tunnel ro- 
tary balance test data (ref. 12). In free spinning 
model tests, prospin controls are employed. Since 
the deflection of control surfaces is not considered in 
the present analysis, the spin-tunnel rotary balance 
test data included in figures 35 to 39 are also taken 
for zero control surface deflections. This kind of com- 
parison is not a true indication of the real situation 
but still is a good representative of the aerodynamic 
parameters dictating the spin modes. 

From figures 45, 47, and 49, which correspond 
to steep or moderately steep spin modes for tails 3, 
4, and 5, respectively, it is observed that for most 
of the coefficients there is a fair agreement between 
the present theory and spin-tunnel rotary balance 
test data (ref. 12), around the values of oo where the 
steady-state spin modes occur. However, for fiat spin 
modes where the reduced spin rates are high (figs. 46 


and 48), the discrepancies in the estimation of Cy 
and C n are large and should be of concern. There- 
fore, the present theory needs further development, 
particularly with respect to Cy and C n for satisfac- 
tory prediction of flat spin modes. 

Concluding Remarks 

The aerodynamics of a spinning airplane is com- 
plex and is dominated by extensive regions of sepa- 
rated flow fields, mutual interferences between aero- 
dynamic and control surfaces, and the complex effect 
of spin rate. As a result, the aerodynamic coeffi- 
cients of a spinning airplane display a high degree 
of nonlinearity with angle of attack and spin rate; 
this makes the conventional approach based on sta- 
bility and control derivatives invalid for the problem 
of the spinning airplane. This complex problem is 
still not amenable to computational fluid dynamic 
methods. The present theory which is semiempirical 
in nature and is developed on the basis of synthesis of 
the available experimental data from various sources 
is an attempt to fill the void. In this way, this work 
represents a first, unified approach to the problem of 
estimating the aerodynamic characteristics of a light 
spinning airplane. In the course of this investigation, 
it was realized that several lacunae exist in literature, 
and this problem, which is of vital importance to gen- 
eral aviation, needs immediate attention to generate 
an experimental data base, which can give more in- 
sights and help in the development of a theoretical 
model. In view of such a situation, the present theory 
had to be based on many speculations and heuristic 
assumptions to put all the pieces together. The re- 
sults and experience of this study are summarized as 
follows along with the recommendations for future 
work. 

The strip theory approach is capable of predicting 
the aerodynamic characteristics of a steadily spinning 
airplane. With appropriate input data to account for 
the effects of angular velocity in spin, its accuracy 
is significantly enhanced. The present semiempirical 
theory is found to be capable of providing reason- 
able estimates of the aerodynamic coefficients of a 
spinning airplane for steep and moderately steep spin 
conditions. However, the theory needs further devel- 
opment for successful application to flat spin prob- 
lems. 

At high angles of attack a, the rotary balance test 
data indicate that the spin rate has a strong influ- 
ence on normal-force, pitching-moment, and rolling- 
moment coefficients of a spinning airplane. In com- 
parison, the predictions based on strip theory with 
static wind-tunnel data as input differ significantly 
from these test data. The physical mechanism re- 
sponsible for this phenomenon is not clear at this 
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stage as there is no experimental data base which 
records flow pattern or pressure distributions over a 
spinning wing and throws some light on the physics 
of the flow. 

No validated aerodynamic models for stalled spin- 
ning wings exist. McCormick (NASA CR-165680) 
hypothesizes incremental changes in the aerodynamic 
forces and moments caused by the spanwise pressure 
gradient on the lee side of the wing due to centrifu- 
gal force. Although the inclusion of the increments 
improves the comparison of strip theory with exper- 
imental test data, there is a dearth of evidence nec- 
essary to validate the underlying assumptions. Lo- 
cal flow measurements of stalled spinning wings are 
needed to enable the formulation of a comprehensive 
theory. 

A semiempirical method is developed for the pre- 
diction of body aerodynamic characteristics at high 
angles of attack and sideslip. This approach gives 
the required static input data in the strip theory for 
the spinning airplane. 

All the mutual interference effects between vari- 
ous components are not considered except the shield- 
ing of vertical tail surface by the wake of horizontal 
tail. A simple procedure for estimating the shielded 
area of the vertical tail is presented. Increments to 
aerodynamic coefficients on account of the secondary 
flow within the shielded part of the vertical tail are 
estimated by a simple analysis of wake flow impinge- 
ment on the vertical tail. 

Some of the mutual interference effects between 
body, wing, and tail surfaces of a spinning airplane, 
which are not considered herein, are found to be large 
and significant. Of particular concern is the interfer- 
ence effect between aft body and horizontal tail sur- 
faces, which can have a strong influence on side force 
and yawing moments. Further studies are necessary 
in this direction to gain a better understanding of 
these complex effects. 

It is recommended that further research be initi- 
ated to study the aerodynamics of the spinning air- 
plane at high angles of attack and spin rate. Flow 
visualization and pressure distribution measurements 
should be carried out along with direct force and 
moment coefficients so as to increase the process of 
understanding and development of a comprehensive 
theory for a spinning light airplane. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
November 5, 1987 
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Table I. Dimensional Characteristics of Full-Scale Model A Airplane 

Wing: 

Span, m 7.46 

Area, in 2 9.11 

Root chord, m 1.22 

Tip chord, m 1.22 

Mean aerodynamic chord, m 1.22 

Aspect ratio 6.10 

Dihedral, deg 5.0 

Incidence at — 

Root, deg 5.0 

Tip, deg 5.0 

Airfoil section Modified NACA 642-415 

Fuselage: 

Overall length, m 5.39 

Distance from leading edge to airplane center of gravity, m 1.91 

Maximum width, m 1.09 

Maximum height, m 1.15 

Maximum cross-sectional area, m 2 1.16 

Volume, m 3 3.48 

Horizontal tail: 

Span, m 2.34 

Incidence, deg 3.0 

Root chord, m 1.10 

Tip chord, m 0.51 

Mean aerodynamic chord, m 0.84 

Aspect ratio 2.91 

Tail length (distance of center of gravity to 0.25 ch), In, for tail 3, m 3.48 

Airfoil section NACA 65i-012 

Vertical tail: 

Span, m 1.25 

Root chord, m 1.10 

Tip chord, m 0.51 

Mean aerodynamic chord, m 0.84 

Aspect ratio 1.55 

Tail length (distance from center of gravity to 0.25cy), /y, for tail 3, m 3.48 

Airfoil section NACA 65i-012 
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Table V . Rolling-Moment Coefficients of Airplane Components 
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Figure 1. Schematic diagram of steadily rotating airplane model in spin tunnel. 
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Figure 3. Variation of section lift and drag coefficients on rotatin 
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No. 625 


propeller. (From Himmelskamp, ref. 9.) 






Figure 5. Comparison of strip theory calculations for normal-force coefficient with rotary balance test data. 
(From ref. 13.) 



Figure 6. Comparison of McCormick’s strip theory calculations for normal-force coefficient incorporating radial 
pressure gradient term (ref. 13) with rotary balance test data (ref. 11). 
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01 


Figure 7. Comparison of calculated rolling-moment coefficient based on strip theory and radial pressure 
gradient correction (ref. 13) with rotary balance test data (ref. 11). 
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Right wing 



Figure 11. Schematic diagram showing shielding of vertical tail in right spin. 











Q 


<? 



Side force due to secondary flow effect 


Strip theory 

Strip theory + secondary flow effect 
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Figure 14. Schematic illustration of secondary flow effects over vertical tail in right spin. 
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(a) Normal-force coefficient (from ref. 10). 
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(b) Axial-force coefficient (from ref. 10). 
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(c) Center of pressure (from ref. 33). 

Figure 18. Input static aerodynamic data for wing and horizontal tail. Nr,, = 0.288 x 10 6 . 
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Figure 23. Correlation of side-force data (ref. 10). Nft e = 0.288 x 10 6 . 
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Figure 29. 


Static longitudinal aerodynamic coefficients of body of model A airplane. Nr € — 0.288 x 


10 6 . 
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(c) a = 60°. 


(d) a = 90°. 


Figure 40. Concluded 
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Spin-tunnel rotary balance test data {ref. 11) 
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(a) a = 30°. 


(b) a = 45°. 


Figure 44. Rotary aerodynamic characteristics of BWHV configuration with tail 5. 
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Figure 45. Rotary aerodynamic characteristics of 
model A airplane with tail 3 at moderately flat 
spin attitude, a = 50°. 


Figure 46. Rotary aerodynamic characteristics 
of model A airplane with tail 3 at flat spin 
attitude, a = 80°. 
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Figure 47. Rotary aerodynamic characteristics of 
model A airplane with tail 4 at steep spin 
attitude, a = 35°. 
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Figure 48. Rotary aerodynamic characteristics 
of model A airplane with tail 4 at flat spin 
attitude, a = 80°. 
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